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Abstract 

Recently, Kalkreuter obtained the complete Dirac spectrum for an SU (2) lattice gauge 
theory. We performed a statistical analysis of his data and found that the eigenvalue cor- 
relations can be described by the Gaussian Symplectic Ensemble. Long range fluctuations 
are strongly suppressed: the variance of a sequence of levels containing n eigenvalues on 
average is given by (log n + const.). Our findings are in agreement with the anti-unitary 
symmetry of the lattice Dirac operator for Nc = 2 with staggered fermions. For Nc = 3 
we predict that the eigenvalue correlations are given by the Gaussian Unitary Ensemble. 
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The QCD Dirac operator is of fundamental importance for the calculation of the 
physical properties of QCD. Knowledge of its eigenvalues and eigenfunctions determines 
the propagator, a necessary ingredient for the calculation of hadronic correlation functions. 
In this letter we will focus on the eigenvalues of the Dirac operator and isolate some 
universal properties that can be understood from its symmetries only. One such example 
is already known: the eigenvalues near zero satisfy sum-rules with a generating function 
given by random matrix theories with the global symmetries of QCD 0. This led us to the 
conjecture that the fluctuations of the eigenvalues no more than a few level spacings away 
from zero, over the ensemble of gauge field configurations are given by universal functions 
that can be obtained from a much simpler random matrix theory as well. This raises 
the question whether the eigenvalues in the bulk of the spectrum show such universal 
characteristics as well. 

Recently, in a ground breaking work, Kalkreuter |Q] succeeded to compute the com- 
plete spectrum of the Dirac operator on a reasonable large lattice (12^). His results were 
in complete agreement with an analytical sum rule adding greatly to our confidence in the 
accuracy of his results. Long level sequences have been analyzed before in atomic and nu- 
clear physics and for systems with only a few degrees of freedom 0, ^ . Generically, it was 
found that if the system is classically chaotic, the correlations between eigenvalues with 
the same exact quantum numbers are given by the gaussian random matrix ensembles. 

In this letter we will perform a statistical analysis of the lattice spectra using such 
methods 0. We start from the assumption that the average eigenvalue density, p(A), can 
be separated from the fluctuations of the eigenvalues about their average position. This 
allows us to unfold the spectrum. This is a procedure in which the eigenvalues are rescaled 
according to the average local level density. Formally, the unfolded spectrum {A'„}, with 
average eigenvalue density p'(A) = 1, is given by 



From the analysis of the spectra of the Hamiltonians of classically chaotic systems, we 
have learned that depending on the time reversal symmetry of the system the level corre- 
lations fall into three different universality classes. We want to stress that an anti-unitary 
symmetry determines whether the matrix elements are complex, real or quaternion real ^ . 
The corresponding invariant random matrix ensembles are called the Gaussian Unitary 




(1) 
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Ensemble (GUE), the Gaussian Orthogonal Ensemble (GOE) and the Gaussian Symplec- 
tic Ensemble (GSE), respectively. According to general universality arguments the 
correlations between eigenvalues in the bulk of the spectrum, as opposed to those near 
the edge of the spectrum, are not sensitive to many other details of the random matrix 
model. For example, the eigenvalue fluctuations in the bulk of the spectrum of random 
matrix models 0, with the chiral symmetry of the Dirac operator built in, will be given 
by one of the invariant random matrix models (GUE, GOE or GSE). 

Let us analyze the anti-unitary symmetries of the Euclidean Dirac operator, T> = 
"ilfidn + jfj.An, for fundamental fermions in an arbitrary SU{Nc) background gauge field 
Afj^. For three or more colors there are no anti-unitary symmetries and the matrix elements 
of the Dirac operator T> are complex. For two colors the situation is different. In the 
continuum theory we have 

Crsi^] = 0, (2) 

where C is the charge conjugation matrix (C = 7274), T2 is one of the Pauli spin matrices, 
and K is the charge conjugation operator. Naive lattice fermions also obey this symmetry, 
but for Wilson fermions it is violated by the r— term. For staggered lattice fermions the 
only remnant of the 7— matrices is a phase factor ±1 and instead of (2) we have p|] 

The important difference is that 

{CT2Kf = l, but {T2Kf = -l. (4) 

From a similar analysis of the time reversal operator in quantum mechanics (see 0) we 
conclude that in the continuum theory the matrix elements of the Dirac operator can be 
chosen real whereas for staggered fermions they can be organized into real quaternions . 
Therefore, we expect that the level correlations of the eigenvalues of X>^"^^- are described 
by the GOE, whereas for V^^^^- they are given by the GSE. A necessary condition in both 
cases is that the gauge potential is 'sufficiently random'. For Wilson fermions the matrix 
elements are complex even for two colors, but at present it is not clear whether breaking 
of the anti-unitary symmetry by an irrelevant operator leads to different level statistics. 
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In his work, Kalkreuter gives results for only a few SU{2) gauge field configurations, 
namely for (3 = A/g"^ = 1.8 and (3 = 2.8, both with periodic and anti-periodic boundary 
conditions in the Euclidean time direction (in the spacial directions he uses periodic 
boundary conditions). His configurations have been obtained for dynamical staggered 
fermions with a mass of 0.2. For (3 = 1.8 the theory is certainly in the strong coupling 
phase, but it is generally believed that asymptotic scaling has been reached for (3 = 2.8. 

With the eigenvalues of only one configuration we are unable to perform an ensemble 
average of their correlations. Instead we will perform a spectral average, which, at least 



in random matrix theory can be shown to be equal to the ensemble average |T0[ . For the 
present data we have verified numerically that the eigenvalue correlations do not change 
over the range of the spectrum. 

The integrated level density (see Fig. 1) 

N{X)= f'p{X')dX', (5) 





of the eigenvalues of V^^^^- follows immediately from the 10368 eigenvalues calculated by 
Kalkreuter ||^. Except for a possible nonanalytical behaviour at the ends of the spectrum, 
it is extremely smooth and almost linear. To unfold the spectrum we cut it in pieces of 
500 eigenvalues and fit them by a second order polynomial. We have checked that this 
results in an average eigenvalue density equal to one, and that our results are insensitive 
to the details of the cuts. 

To define our statistics we introduce the quantity 

N{x,x + n) = r^"" p'{X)dX, (6) 

Jx 

where p'(A) is the unfolded spectral density. So, N{x, x + n) is the number of eigenvalues 
in a sequence of length n starting at x. The number variance is defined by 

S2(n) = - E(iV(x„ x, + n)- n)\ (7) 
P i=i 

where the points Xi are regularly spaced such that the sequences [xi,Xi + n] cover the 
complete spectrum. It can be related to the two-point level correlation function which is 
known analytically in random matrix theory (see lll|). For a random sequence of levels 
(Poisson spectrum) it can be shown that Ti2{n) = n (see P]). 



Our results for S2(n) are shown in the upper figure of Fig. 2 (/3 = 1.8) and Fig. 3 (/3 = 
2.8). (To ehminate inaccuracies in the unfolding procedure we excluded 500 eigenvalues 
at the ends of the spectrum.) Long range fluctuations are almost completely absent. 
Instead of a variance of 100 for an average sequence of length 100 we flnd a variance of 
only 0.4, showing the presence of very strong correlations between the eigenvalues. From 
what we have said before, we expect that they can be described by the GSE. Indeed, the 
theoretical result (full curve), for r > 1 given by 

S2(r) = -i^(log(47rr) + ^ + 1 + !^) + O(^), (8) 

shows a perfect agreement with the lattice data. For comparison, we have also given the 
result for the GUE (dashed curve). The result for the GOE is for the most part outside 
the range of the flgure. 

A much smoother statistic is the As-statistic originally introduced by Dyson and 
Mehta [|l^]. It is related to the number variance by ||13| 



Mr) = -J (r' - 2r^s + /)S2(s)cis, (9) 
JO 

and does not receive contributions from the quadratic term in S2(s) so that small unfolding 
errors are eliminated. Asymptotically, for r > 1, one flnds for the GSE 

A3(r) = \^,{r) - (10) 

Results for A3 are shown in the middle flgure of Fig. 2 and Fig. 3. Numerical results 
obtained from the data of Kalkreuter are given by full dots whereas the full curve refers 
to the GSE result and the dotted curve represents the A3 statistic of the GUE. Also in 
this case the curve for the GOE lies well above the GUE. 

Finally, in the lower flgure of Fig. 2 and Fig. 3, we show the distribution P{S) of 
the nearest neighbor spacings S = — X'i\ of the unfolded spectrum. The histograms 
obtained from the lattice data are represented by points. The full lines show the GSE 
result for the nearest neighbor spacing distribution and the dashed curve give the result 
for the GUE. The exact expression is very well approximated by the Wigner surmise which 
for the GSE is given by 

PiS)^^S'e.,{-'^J^). (U) 



For comparison we also show the result for the GUE (dashed curve) which is quadratic 
in S for small S. Concerning the spacing distribution of the GOE we only remark that it 
starts out linearly and deviates even more from the GSE. The nearest neighbor spacing 
distribution contains information on short range correlations in the spectrum. We hope 
that we have convinced the reader that they are also given by the GSE. 

We have also analyzed the spectra that Kalkreuter obtained for periodic boundary 
conditions in all directions and found identical results for the eigenvalue correlations. 

In conclusion, we have found that the fluctuations of the eigenvalues of the staggered 
Euclidean Dirac operator for SU{2)-coIot can be described by the Gaussian Symplectic 
Ensemble both for (3 = 1.8 and (3 = 2.8. For the continuum theory the anti-unitary 
symmetry is different and we expect fluctuations according to the Gaussian Orthogonal 
Ensemble. In ||2|, we have shown that anti-unitary symmetries are essential in de- 
termining the Goldstone sector of the theory. Because the massless sector for staggered 
fermions differs from the continuum theory, it is not surprising to find level correlations 
that belong to a different universality class. It would be very interesting to analyze the 
fate of the anti-unitary symmetries and eigenvalue correlations in the continuum limit of 
the lattice gauge theory. Lattice simulations with larger values of (3 and bigger lattices 
are required to investigate this point. 

For three or more colors the anti-unitary symmetries are broken both in the continuum 
theory and for staggered lattice fermions. We predict that in this case the eigenvalue 
correlations are given by the GUE. It would be worthwhile to obtain the complete Dirac 
spectrum also for this case. 

A final point of interest we want to mention is the fate of level correlations during the 
chiral phase transition. From solid state physics we know a delocalization transition 
is associated with a transition in the level statistics which raises the hope that such 
phenomena can be seen in QCD as well. 
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Figure Captions 



Fig. 1. The integrated eigenvalue density A'"(A) of the Euchdean lattice Dirac operator for 
SU (2)-color with staggered fermions. Each curve is the result for one equilibrated gauge 
field configuration with (3 as indicated in the figure. 

Fig. 2. Results for the level statistics T,2{n) (upper), A3(n) (middle) and the nearest 
neighbor spacing distribution P{S) (lower) of the eigenvalues shown in Fig. 1 for /3 = 1.8. 
The random matrix results for the GSE and the GOE are represented by the full and the 
dashed line, respectively. 

Fig. 3. Results for the level statistics ^2(71) (upper), A3(n) (middle) and the nearest 
neighbor spacing distribution P{S) (lower) for /3 = 2.8. For further explanation see the 
caption of Fig. 2. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9501025vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9501025vl 



This figure "figl-3.png" is available in "png" format from: 
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